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Under a two-year program at the General Electric Company, a series of three- 
dimensional inelastic structural analysis computer codes were developed and 
delivered to NASA.* A summary of the present state of the capabilities of these 
computer codes is presented here. 

The objective of this program was to develop analytical methods capable of 
evaluating the cyclic time-dependent inelasticity which occurs in hot section engine 
components. Because of the large excursions in temperature associated with hot sec- 
tion engine components, the techniques developed must be able to accommodate large 
variations in material behavior including plasticity and creep (ref. 1). To meet 
this objective. General Electric developed a matrix consisting of three consti- 
tutive models and three element formulations. A separate program for each combina- 
tion of constitutive mode - element model was written, making a total of nine 
programs. The source codes of the nine programs range in size from 7300 lines for 
the Bodner/twenty node to 19 000 lines for the Haisler and Allen/nine node. Table I 
shows the length of each source code. All of the codes were given a stand-alone 
capability of performing cyclic nonlinear analysis. 

The three constitutive models consist of a simple model, a classical model, and 
a unified model. In an inelastic analysis, the simple model uses a bilinear stress- 
strain curve to determine the plastic strain and a power law equation to obtain the 
creep strain. The second model is the classical model of Haisler and Allen (ref. 2). 
The third model is the unified model of Bodner, Partom, and Partom (ref. 3). The 
attributes of the three constitutive models are listed in Table II. All of the 
models were programmed for a linear variation of loads and temperatures with the 
material properties being temperature dependent. 

The three element formulations used are an eight-node isoparametric shell 
element, a nine-node shell element, and a twenty-node isoparametric solid element. 

The eight-node element uses serendipity shape functions for interpolation and 
Gaussian quadrature for numerical integration. Lagrange shape functions are used in 
the nine-node element. For numerical integration, the nine-node element uses 


1-Work done under NASA contract NAS3-23698. 




rX&VcD 


267 



Simpson's rule. The twenty-node solid element uses Gaussian quadrature for 
integration. The other attributes of these elements are listed in Table III. 

For the linear analysis of structures, the nine codes use a blocked-column 
skyline, out-of-core equation solver. To analyze structures with nonlinear material 
behavior, the codes use an initial stress iterative scheme. Aitken's acceleration 
scheme was incorporated into the codes to increase the convergence rate of the 
iteration scheme. 

The ability to model piecewise linear load histories was written into the 
codes. Since the inelastic strain rate can change dramatically during a linear load 
history, a dynamic time-incrementing procedure was included. The maximum inelastic 
strain increment, maximum stress increment, and the maximum rate of change of the 
inelastic strain rate are the criteria that control the size of the time step. The 
minimum time step calculated from the three criteria is the value that is used. 

In dynamic analysis, the eigenvectors and eigenvalues can be extracted using 
either the determinant search technique or the subspace iteration method. These 
methods are only included with those finite-element codes containing the eight-node 
shell element. 

The nine codes have been compiled on the CRAY-1 machine at NASA Lewis Research 
Center. Table IV shows those features contained in the codes that have been 
checked. Additional work has to be done in examining the other features contained 
in these codes. Also, test problems need to be analyzed so that the performance of 
these codes can be examined. 
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TABLE I. - APPROXIMATE LENGTH OF THE NINE 
SOURCE CODES 


Constitutive 
model s 

Element formations 

20-Node 

8-Node 

9-Node 

Length of code, number of lines 

Simpl e 

8300 

13 800 

17 900 

Classical 

9200 

16 300 

19 000 

Uni f i ed 

7300 

13 800 

17 600 


TABLE II. - CONSTITUTIVE MODELS 


Simple model 

Classical model 

Unified model 

Uncoupled plasticity 
and creep 

Uncoupled plasticity 
and creep 


Plasti ci ty 

PI asti ci ty 

No yield surface 

Isotropic hardening 

Combined isotropic and 
kinematic hardening 

Isotropic hardening 

Piecewi se 1 i near 

stress-strain curves 

Piecewise linear 

stress-strain curves 


Prandtl-Reuss 
flow rule 

Modified Prandtl-Reuss 
flow rule 



Noni sothermal 

Noni sothermal 

Noni sothermal 

Creep 

Steady state 
Isotropic hardening 
Prandtl-Reuss 
flow rule 
Noni sothermal 

Creep 

Steady state 
Isotropic hardening 
Prandtl-Reuss 
flow rule 
Noni sothermal 

Second-order Adams- 
Moulton rule used 
to carry out 
i ntegrati on 
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TABLE III. - ELEMENT FORMULATION 


8-Node shell 

9-Node shell 

20-Node solid 

Five degrees of freedom 
3 Displacements 
2 Rotations 

Five degrees of freedom 
3 Displacements 
2 Rotations 

Three degrees of 
freedom 

3 Displacements 

Serendipity shape 
functions 

LaGrange shape functions 


No rotational stiffness 
about the normal to the 
mid-surface; deleted 
prior to assembly 

Rotation about the normal 
to the mid-surface is 
treated as a prescribed 
di spl acement 


Isotropic or orthotropic 
elastic properties 

Isotropic or orthotropic 
elastic properties 

Isotropic or 
orthotropic elastic 
properties 

Surface, line, nodal, 
RPM, thermal , and 
gravity loads 

Surface, line, nodal , 
RPM, thermal , and 
gravity loads 

Surface, nodal, 
RPM, thermal , and 
acceleration loads 

Prescri bed 
di splacements 

Prescribed displacements 

Prescri bed 
displacements 

Gaussian quadrature 
used for numerical 
integration 

Simpson’s rule used for 
numerical integration 

Gaussian quadrature 
used for numerical 
integration 
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TABLE IV. - CODE FEATURES 3 


Feature 

Simple model 

Classical model 

Unified model 


8-Node 

9-Node 

20-Node 

8-Node 

9-Node 

20-Node 

8-Node 

9-Node 

20-Node 

Free format data input 
Global coordinate system 

X a 

X 


(■ 

X 

X 

X 

X 

X 

Cartesian 

X 

X 



X 

X 

X 

X 

X 

Spherical® 

— 

— 


1 

«... 









Cyl indri cal® 

Local coordinate system 

— 1 

— 

Bill 

— 

— 

1 

— 

— 

— 

Cartesian 

X 

X 

X 

X 

X 

X 

X 

X 

... 

Spherical® 

— 

— 

— 

— 

— 

— 







Cyl i ndri cal ® 

— 

— 

— 

— 

— 

— 







Automatic generation of 

— 

■ 

N/A c 

— 

. 

N/A 





N/A 

nodal coordinates 
Automatic generation of 

X 

X 

N/A 

X 

X 

N/A 

X 

X 

N/A 

element connectivities 
Restart capability 

_ 








B 

Dynamic allocation 

X 

X 


X 

X 

■ 

X 

X 


Blocked column skyline 

X 

X 


X 

X 


X 

X 


equation solver 
Initial stress iterative 



1 






■9 

scheme 

Aitken's acceleration 



X 



X 



X 

scheme 










Oynamic time 
incrementing 
Convergence criteria 



X 



X 



X 

Effective plastic 

— 

— 

X 

— 

— 

X 

— 

— 

X 

strai n 

Effective stress 





X 


— 

X 



_____ 

X 

Dynamic analysis 
Eigenval ue 

X 

N/A 

N/A 

X 

N/A 

N/A 

X 

N/A 

N/A 

Ei genvector 

X 

N/A 

N/A 

X 

N/A 

N/A 

X 

N/A 

N/A 

Material change option 
Element removal option 


IBb 

— 

:: 

“ — 

— 

Bj 

B 

B 

Loads 



X 

X 

X 

X 




Midsize node generation 




— 








Skewed coordinate system 



X 

X 

X 

X 

■ 

life ' 


Orthotropic orientation 



X 

X 

X 

X 




def i ni ti on 

■■ 

| 





Hi, 

■■ 

■■ 


3 X = Feature has been checked. 

®Cylindrical and spherical coordinate systems have not been checked. 
C N/A = Not a feature of this particular code. 
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